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Abstract 

In the present paper, a class of stochastic Runge-Kutta methods containing the se- 
cond order stochastic Runge-Kutta scheme due to E. Platen for the weak approxi- 
mation of Ito stochastic differential equation systems with a multi-dimensional 
Wiener process is considered. Order one and order two conditions for the coef- 
ficients of explicit stochastic Runge-Kutta methods are solved and the solution 
space of the possible coefficients is analyzed. A full classification of the coefficients 
for such stochastic Runge-Kutta schemes of order one and two with minimal stage 
numbers is calculated. Further, within the considered class of stochastic Runge- 
Kutta schemes coefficients for optimal schemes in the sense that additionally some 
higher order conditions are fulfilled are presented. 
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1 Introduction 



Recently, the development of numerical schemes for strong as well as weak 
approximation of stochastic differential equations (SDEs) has focused amongst 
others on Runge-Kutta type schemes pl|2 5 6 7 8 }9|10|ll|ll2|ll3j . This is due 
to the increasing complexity of stochastic Taylor expansions and the desire 
to avoid derivatives in higher order approximation schemes. In section El a 
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class of stochastic Runge-Kutta (SRK) methods due to Rofiler [TC)|11|I12] for 
the weak approximation of Ito SDE systems with a multi-dimensional Wiener 
process is considered. This class contains as a special case the second order 
SRK scheme proposed by Platen [5] as well as the class of SRK methods pro- 
posed by Tocino and Vigo-Aguiar [13] . Order conditions for coefficients of the 
SRK methods have been calculated by applying the multi-colored rooted tree 
analysis due to Rofiler [T0fll|ll2] . In contrast to earlier work on this topic, the 
aim of the present paper is to analyze these order conditions with the objec- 
tive to determine a full classification of the coefficients for this class of SRK 
methods. A full classification for order one SRK schemes with s = 1 stage 
and for order one SRK schemes with deterministic order two for s = 2 stages 
as well as for second order SRK schemes with s = 3 stages is calculated in 
section [31 Further, some optimal schemes are derived from this classification 
in section H] by taking into account additional higher order conditions. Their 
performance is studied by some numerical examples in section 

We denote by (X t )t&i the solution of the rf-dimensional Ito SDE defined by 

dX t = a(t, X t ) dt + b(t, X t ) dW t , X t0 = x , (1) 

with an m-dimensional Wiener process (W t )t>o and I = [to,T]. We assume 
that the Borel-measurable coefficients a: IxR d ^R d and b : I x R d — > R dxm 
satisfy a Lipschitz and a linear growth condition such that the Existence and 
Uniqueness Theorem [5] applies. In the following, let V(t, x) = (6 lJ (t, x))\<i<d G 
R d denote the jth column of the diffusion matrix b(t, x) for j = 1, . . . , m. 

Let a discretization I h = {t , ti, . . . , t]y} with t < t\ < . . . < tjq = T 
of the time interval I = [to,T] with step sizes h n = t n +\ — t n for n = 
0, 1, . . . , AT — 1 be given. Further, define h = maxo< n <7v h n as the maximum 
step size. Let C l P (R d , R) denote the space of all g E C l (R d , K) fulfilling a poly- 
nomial growth condition and let g G Cp l (I x ]R d ,M) if g(-,x) G C k (I, R) and 
g(t, •) G C l P (R d , R) for all t G / and x G R d 0. 



Definition 1.1 A time discrete approximation Y = (Y t ) te j h converges weakly 
with order p to X as h — > at time t G Ih if for each f G Cp (M. d , R) exist 
a constant Cf and a finite 5q > such that 

\E(f(X t ))-E(f(Y t ))\<C f h? (2) 

holds for each h G ]0, 5q[ . 
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2 Stochastic Runge— Kutta Methods 



We consider stochastic Runge-Kutta methods as proposed in [T0|llfl2] for the 
weak approximation of SDE (JTJ). Therefore, the d-dimensional approximation 
process Y of an explicit s-stage SRK method is defined by Y to = xq and 

Yt n+1 = Y tn + £ a, a(tn + cf ] h n , HP) h n 

i=l 
s m 

i=l k=l 
s m 

+ ££/3f ) ^n + C?K,tff ) )^ 

i=l k=l {o) 
s m 

+ £ £ A (3) b k (t n + cf ] h n , hP) i {k) 

i=l k,l=l 
k^l 
s m 

+ EEA (4) ^ + c«\f)^ 

i=l k,l=l 
k^l 

for n — 0, 1, . . . , N — 1 with stage values 

H? ] = Y tn + £ 4? a(t n + cf hn, Hf) h n 

3=1 
i—l m 

+ Y,Y,B [ S ) v\t n + cfh n ,Hf ) )i {T) 

j=l r =l 

H\ k) = Y tn + l £ Af) a(t n + cfh n , Hf) h n 

3=1 

+£5^ b k {t n + c fh n ,Hf ] )jh n 

3=1 

H\ k) = Y tn + £ 4? <tn + cfh n , Hf) h n 

3=1 

+ £4° b\t n + cfh n ,Hf ) )Jh n 

3=1 

for i = 1, . . . , s and k = 1, . . . , to. Here, a, . . . , /3 (4) , c (g) G R s and A^, 
5(a) e k sxs for < q < 2 with A§ } = B[f = for j > i are the vectors 
and matrices of coefficients of the SRK method. We choose = A^ q 'e for 
< q < 2 with a vector e = (1, . . . , 1) T [TO] . In the following, the product 
of column vectors is defined component-wise. The coefficients of the SRK 
method (02) are determined by the following Butcher tableau: 
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c (0) 




5(0) 




c (l) 




5(1) 




c (2) 


^(2) 


£(2) 






T 

a 










/3( 3 ) r 


/3W T 



The random variables of the SRK method are defined by three-point dis- 
tributed random variables with P(I[ r ) = ±^3 h n ) = | and P(/( r ) = 0) = |. 
Further, I/ k) ^ = \{I(k)I(i) + V k ' 1 ). The V k ' 1 are independent two-point dis- 
tributed random variables with P(V k ' 1 = ±h n ) = | for I = l,...,k — 1, 
yk,k _ _y l ^ anc j yk,l _ _yl,k for Z = A; + 1, . . . , m and fc=l,...,m|. 

By the application of the multi-colored rooted tree analysis [T0~fT2"] . order 
conditions for the coefficients of the SRK method ([3]) can be easily deter- 
mined. As a result of this, the following Theorem 12.11 due to Rofiler [12] gives 
order conditions for the SRK method (J3J) up to order two. 

Theorem 2.1 Let a\b ij e C 2 /(I x R d ,R) for 1 < i < d, 1 < j < m. If the 
coefficients of the SRK method (TJ|) fulfill the equations 

1. a T e = 1 2. /3 (4)T e = 3. /3 (3)T e = 

4. (/3« T e) 2 = l 5. ^ T e = 6. 0^ T B^e=Q 

7. p® T B®e = 

then the SRK method converges with order 1 in the weak sense. In addition, 
if a\ b ij E C 3 /(I x R d , R) for 1 < i < d, l<j<m and if the equations 



8. 


a T A^e = \ 


9. 


a T (5(°' e ) 2 = | 


10. 


{^ T e){a T B^e) = \ 


11. 


(^ T e)((3^ T A^e) = 


12. 


^) T A m e = o 


13. 


pw T B m e = i 


14. 


/3 W T jB (2) e = x 


15. 


(/3« T e)(/3« T (sW e ) 2 ) 


16. 


(^ T e)(^ T (B^e) 2 ) = \ 


17. 


^(i) T ( jB (i)( jB (i) e )) = 


18. 


^) T (5( J) (B (1) e)) = 


19. 


/3( 3 ) T (A( 2 )(s(°) e )) = 


20. 


£(i) T (A (D (jB (o) e)) = 


21. 


a T (jB (0 )(jB (l) e))=0 


22. 


pW T A^e = 


23. 


^(4) T A (2) e = o 


24. 


(3 {1)T ((A^e)(B^e)) = 


25. 


/^((A^eX^e)) = 


26. 


^(4) T (A (2) (jB (o) e)) = 


27. 


/3(2) T ( A (i) (jB (o) e)) = 
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ZO. 


/ o(2) T M (l)/ R (0) n2n _ n 


9Q 


/ o(4) T M (2)/ R (0) n2n _ n 


ou. 


tf(3) T '( D(2) ( 4(1) \\ _ n 




bW t (bW(a(V p )) - 


32. 


/ 3(2) T ( jB (D e )2 = 


33. 


^(4) T (S (2) e) 2 = q 


34. 


/3 (4) r ^(2)^(l) ^ = g 


35. 


R^ T (B^(B^e)) = 


36. 


^(l) T ( S W e )3 = 


37. 


^(3) T ( S (2) e )3 = 


oo. 


rW 7 I rW f rW ^^2n _ n 
p l^ij (^ij ej j — u 


oy. 


o(3) T ( R (2) / R (l) \2\ _ n 




40. « T ((5<°) 


e)(B^(B 


We))) = 



41. (3^ T ((A ( - 1 \B^e))(B^e))=0 

42. /^((^(E^X^e)) =0 

43. /3( 1 ) T (A( 1 )( J B(°)( jB (i) e ))) = 

44. ^ T (A^(B^(B^e))) = 

45. /3« T (S«(A«(S(°)e))) = 

46. l3^ T (B^(A^(B^e))) = 

47. /3^ T ((B^e)(B^(B^e))) = o 

48. ^ T ((B^e)(B^(B^e)))=0 

49. /3W T (sW(BW(BWe))) = 

50. /3^V (2) (£ (1) 0B (1) e))) = O 

are fulfilled then the stochastic Runge-Kutta method |3j) converges with order 2 
in the weak sense. 



In the case of m > 1 one has to solve 50 non-linear equations in order to 
calculate coefficients for an order two SRK method ([3]). However, in the case 
of m = 1 these conditions are reduced to 28 equations which have to be 
solved [TT] . Thus, the analysis of the space of all admissible coefficients is not 
an easy task. It turns out that explicit order one SRK methods need at least 
s — 1 stage while order two SRK methods need s > 3 stages. This is due to 
e.g. conditions 6. and 15., which can not be fulfilled in the case of s < 2 stages 
for explicit order two SRK methods. In the following, we distinguish between 
the stochastic and the deterministic order of convergence. Let Ps = P denote 
the order of convergence of the SRK method if it is applied to an SDE and let 
Pd with pr> > ps denote the order of convergence of the SRK method if it is 
applied to a deterministic ordinary differential equation (ODE), i.e., SDE (pQ) 
with 6 = 0. We also write (pd,Ps) m the following [TT]rr2"] . 
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3 Parameter Families for SRK Methods 



3. 1 Coefficients for SRK Methods of Order (1, 1 ) 



First, we analyze explicit SRK methods (J3J) of order pc = p$ = 1 with s = 1 
stage. Considering the order one conditions 1.-7. in Theorem 12. 1\ the corre- 
sponding coefficients are uniquely determined for C\ G {—1,1} by 

«i = l, rf^cx, (3{ 2) = 0, /3f } = 0, /3i 4) = 0. (4) 

The resulting class of SRK schemes coincides with the well-known Euler- 
Maruyama scheme. 



3. 2 Coefficients for SRK Methods of Order (2, 1 ) 



Next, we consider the case of s = 2 stage explicit SRK methods (jSJ). As already 
mentioned in Section [2], it is not possible to attain order ps = 2. However, we 
can find some SRK methods of order pd = 2 and ps — 1 corresponding 
to the following parameter family: From condition 1. of Theorem 12.11 follows 
oil = 1 — a 2 and taking into account the order 2 condition 8. we obtain 
a 2 = —joj f° r ^21^ 7^ 0- Further, condition 2. yields ^ = — /3 2 , condition 

3. results in (3^ = — /3 2 and condition 5. is fulfilled if j3[ 2 ^ = — /3 2 while 
condition 4. holds for t3^ = c\ — /3 2 with c\ G { — 1,1}. Finally, considering 
condition 6. we need that /3 2 = or = and for condition 7. analogously 

(3) (2) 

that /?2 = or 5 21 = hold. Thus, this class of SRK methods is determined 
by 



a 



I- A. JL 

2c 2 2c 2 





c 2 

o' 

c 3 



— C4 C4 
C6 -C 6 

"0 

c 8 

" 

C10 



A (2) 



c 5 — c 5 
C 7 -c 7 

"0 

c 9 

' 

en 



(5) 



for Ci G { — 1, 1} and c 2 , . . . , Cn G M with c 2 7^ 0, c 4 c w = and c 6 Cn = 0. 
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3. 3 Coefficients for SRK Methods of Order (2,2) 



Now, we consider explicit SRK methods (EJ) of order p D = p s = 2 with 

s = 3 stages. Then, the SRK schemes of the class under consideration are 

completely characterized by the following families of coefficients which follow 

from the order conditions in Theorem 12.11 Due to conditions 13. and 32. we 

need 7^ and from conditions 15. and 36. follows /3<p 7^ 0. Thus, there 

exist no SRK schemes of the considered class attaining order = ps = 2 

with less than 3 stages. Now, by condition 17. follows that = and we 

deduce from 6., 15. and 36. that B^i = — -B21 7^ 0. Analyzing the weights, 

(2) 1 (2) 1 
we calculate from conditions 5., 13. and 32. that (3 2 = — m-, P3 = jtj 

2B 21 2B 21 

(2) 

and (3\ = 0. For c\ G { — 1, 1} we obtain from conditions 4., 6. and 15. the 
weights = d - -&J- and & ] = 0^ = -3^. Now, due to 24. and 

11. we need that A$ = (B^) 2 and = (B^) 2 - A$. Applying now 
conditions 3., 7., 16. and 37. we conclude that B$ + B$ + B$ = 0, B$ ^ 0, 
B { £ ^ B { £ + B® and that /?{ 3) = $ ] = and /# = — 

Further, we can now determine the remaining weig hts as /3{ 4) = 0, (5 { 2 ] = -4- 

2B 21 

and /3$ = W from conditions 2., 14. and 33., and we have ot\ = 1 — ct 2 ~ a 3 

2B 21 

(2) 

due to condition 1. Then, we can consider condition 18. which needs B% 2 = 



(2) 
32 

(2) _ o(2) 



and we thus get with the previous considerations that finally £% = — 

1(2: 

l 21 



(2) 

has to be fulfilled. Now, we obtain from conditions 12. and 23. that AX-,' = 



(2) (2) 

and that A\ 2 = —^31 has to be fulfilled. Continuing in this manner, we have 
to distinguish the following cases: 

(A) For = 0, the parameter family is given by a\ = a-i = \ and B 2 } = ci, 
which follows from conditions 1., 9. and 10. Further, we calculate from 
condition 8. that A 2 ± = 1, from 20. that A% 2 = and from condition 19. 



that A£ 2 > = 0. 



|(2) _ ^ 

(B) For «3 7^ 0, condition 21. yields now that B^ 2 = and we have to consider 
the following cases: 

(a) For B 2 \ = it follows from conditions 9. and 10. that B^ = c\ and 
a 3 = 5. Thus, by condition 1. it follows immediately that a% = 5 — «2- 

(i) If A^i = then condition 8. implies that in this case A^ 2 = 
1 - A { S has to be fulfilled. 

(0) 1-A (0) -A (0) 

(ii) If A 2l 7^ then condition 8. yields that a\ — | o 3 * (0) 32 and 



l(0) A (o) 



z/i 21 



«2 = 2 >) 



(b) For B^ 7^ 0, we calculate from condition 20. that A32 = and from 

(2) (2) 
condition 19. that A\ 2 = which implies that also A 31 = due to 
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^32 

and 
CD 



—A$L. Now, with k = a 2 a 3 (2a 2 + 2a 3 
10. we have to consider the following cases: 
In the case of k > 0, a 2 ^ —a?,, «2 / and a 2 ^ ±y/n it follows 



1) and conditions 9. 



that 



202(02+^3) 



-"31 



203(02+03) ' 



If a 2 = and 0:3 = | then we can conclude that B^ 



to hold. 

For a 2 = 



and B 



(o) 

31 



a 3 and a 2 7^ — | it follows that 



(hi) 

Due to condition 8. it follows that A 



ci 



ci has 



'i + 1 } 

^2 ^ 4o 2 ^ 



1 ■ 

4o2 • 



has to be fulfilled. 



(0) 
31 



l-2o 2 A 

2^7 



(0) 



- A 



(o) 

32 • 



Finally, one can easily check that all remaining conditions of Theorem 12.11 
which have not been mentioned explicitly in our analysis, are fulfilled by each 
parameter family and thus do not contribute any further restrictions for the 
coefficients. 



Summarizing our results, we have the following classification for the SRK 
schemes of order p^ = p$ = 2 for the considered class with s = 3 stages: For 
Ci G { — 1, 1} and c 2 , c 3 , c 4 , c 5 G WL with c 3 7^ and c 4 7^ it holds 



A <2) 



Cl 







(2) J 



Cl 






c 2 c 2 






c 5 -C5 



5 (2) 



[° i 


1 

2c 3 . 




1 

2c 4 


" 


0" 




c 3 







-c 3 


0_ 




" 


0" 




c 4 







-c 4 








Now, the following cases are possible: 

In the case ([A]) , we get with c 2 = c 5 = in (JBJ) 



that 



a 



1 1 

2 2 







A (o) 



"0 





0" 




"0 





0" 


1 








, £?<°) = 


Cl 





























with Ai 0) 



4(0) _ R (o) _ R (o) 

*31 — ^32 — n 31 ~ n 32 

for the scheme due to a 3 = 0. 



(6) 
(7) 

(8) 
(9) 



(10) 



because these coefficients are not relevant 



For the case ()Ba.i|) we get with cq, c-j G R the coefficients 



a 







"0 





0" 




"0 





0" 




, A<°> = 











, B<°> = 















.°7 


1-C 7 







c i 









Published in Mathematics and Computers in Simulation 77 (2008) no. 4, pp. 408-420, 



doi: 10.101 6/j.matcom.2007. 04.016 



Considering the case ( IBa.iij) we obtain for cq, c 7 , c$ G E with c% 7^ that 









"0 





0" 




"0 





0" 


T 

a = 


1 1 — C7 — Cg 1 — C7 — C8 1 

2 2c 6 2c 6 2 


, A<°> = 


c 6 








, = 

















_ C 7 


















• (12) 



Next, we have the case (IBb.ij) with c 2 = c 5 = in c 6 , c 7 , c 8 , c g G E and 

with c 6 and c 6 7^ — C7 7^ 0. Then, it holds with k = c 6 c 7 (2c 6 + 2c 7 — 1) and 



A = ±^£8 for c 6 ^ ±y/K and « > that 















0" 




T 

« = 


1 - c 6 - c 7 c 6 c 7 


, A<°> = 


c 8 
















A - Cg 


Cg 










Ci _£6Tv^H_ 

2 c 6 (c 6 +c 7 ) 

C7±y/it 
2 C7(C6+C7) 



0" 



(13) 



The case ( IBb.iil) yields for c$, c 7 , c 8 G E with c 8 7^ and c 2 = C5 = in (JE!)-© 
the coefficients 















0" 




"0 





0" 


T 

a = 


1 i" 

2 u 2 


, A<°> = 


c 6 








, fl<°> = 


c 8 














1-C 7 


c 7 







. c i 









(14) 



Finally, we have the case fIBb.iiij) for c 6 ,c 7 ,c 8 G E with c 6 ^ {— 1,0} and 
c 2 = c 5 = in ©-© which leads to 



a 



1 c 6 -c 6 



A (0) 





C7 



0' 




-2c 6 ~ C 8 C « U . 



B (0) 







a 1 + J- 1 

£1(1 _ _L' 

2 \ 2c 6 < 



0' 



(15) 



3.4 Coefficients for SRK Methods of Order (3,2) 



If we consider the classification of the coefficients for explicit SRK methods, 
we can see from Section 13.31 that in some of the resulting cases there are still 
degrees of freedom in choosing the coefficients for a and A^°\ Therefore, we 
analyze now the classification for explicit SRK methods ([3]) with s = 3 stages 
of order pr, = 3 and ps = 2. Thus, we additionally have to take into account 
the well known deterministic order 3 conditions [3]H] 

a T {A^ef = \, (16) 
a T {A^{A^e)) = -. (17) 
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Clearly, these conditions can not be fulfilled in case (jAj) where a% = as well 
as in case fIBa.ij) due to 4? = and in case fIBb.iij) due to the restrictions for 
a and A^°\ However, in the case of parameter family (IBa.iij) we obtain from 
(HSI) and flUD an SRK method of order (3,2) if in JEJ it holds 



-c 6 ±-y9c| 



36c 6 



24 

3c 6 



c 8 



1 

3c 6 



For the parameter families in case (IBb.ij) and fIBb.iiij) we have to distinguish 



the following three possibilities due to condition 8. of Theorem 12. II and due to 
(USD: 

a ^ nn _ 3 AO) _ 2 AO) _ _A0) 

u\ „ _ 3 n AO) _ 2 AO) _ 2 AO) 

\ _ i 2-3(4°' +A^° 2 ) ) _ 1 34i'-2 - f AO) , .(0) / 

C / tt2 ~~ 6 A (0) (A (0) AO) 7M7> °3 — 6( ,(0) ,«,,,(0) .(0) ,(0K 11 ^31 + A}2 7* 

A 21 ^21 _/1 31 ^^32 -> 1.^31 +/1 32 >( A 2X _/i 31 _/i 32 <> 

4? + and 4? + 4? 1 0. 

If we consider now the case of the parameter family (1Bb.il) then the conditions 
(USD and ([TID are fulfilled for [a} if 

c 6 = l c 7 = i-, c 8 = 2 -, (19) 

with c 7 i {-|,0, §}U] - |,0[ in flU). Further, the conditions (USD and flT7]) 
are also fulfilled in the case (1Bb.il) combined with |b]) if 

3 1 1 2 

Ce = 7 ~ T~ ' C7 = F" ' c § = o ' ( 20 ) 

4 4Cg 4c 9 3 

with c g 7^ and for c 6 £]0,t[u]|,|[ in ([TBI . Finally, the considered parameter 
family (IBb.il) fulfills the conditions (|T6|) and (|T7j) also for the case [cj) if 

1 2-3A 1 3c 8 -2 A(c 8 - A) 

C6 6c 8 (c 8 -A)' ° 7 6A(c 8 -A)' C9 (3c 8 -2)c 8 ^ 

in (USD for A G R with A £ {0, §, c 8 , | - c 8 }, (A - l)c 8 ^ A 2 - §, c 8 £ {0, §} and 
with A < | if c 8 = 1; ^ < A < § if | < c 8 < 1 holds; | < A or A < ^ 

if < c 8 < | holds and with A < § or A > 3 ^~^ if c 8 < or c 8 > 1 holds. 

Note that A = is thus automatically fulfilled in (TT3|) . 



Finally, we consider the case of parameter family ( 1Bb.iii|) which fulfills the 



additional order three conditions ( !T6|) and ( TT7D in case [a| for C6 = f and 
c 7 = | as well as in case[c]) with cq = Ac ^_i if c 7 ^ { — |, 0, |} in (fT5|) . For case 

ED there exists no solution. 
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4 Optimal SRK Schemes 



In the present section, coefficients for the SRK method of different orders 
of convergence are presented. Due to some degrees of freedom in choosing the 
coefficients, we consider additional conditions in order to specify the free pa- 
rameters of the SRK scheme. Clearly, these additional conditions need not ne- 
cessarily be fulfilled for the desired order of convergence. However, coefficients 
fulfilling also higher order conditions yield SRK schemes with the objective to 
obtain smaller error constants and we call them optimal SRK schemes in the 
following. 



4-1 Coefficients for Optimal SRK Schemes of Order (2,1) 

For SRK methods of order pu = 2 and p$ = 1, we need 2 stages for the 
drift part, however only one stage is needed for the diffusion. Therefore, we 
consider only the case of C4 = . . . = en = in Next, we want to specify 
C2 and C3. Therefore, we consider additional order conditions which need not 
necessarily be fulfilled for order (2, 1) schemes. Taking into account condition 
9. yields c 2 = c\. From condition 10. it follows that c 2 = C\ C3. Further, one can 
consider the deterministic order 3 conditions (1TB]) and ffT7]) [3.4], whereas only 
(TIB"]) can be fulfilled which yields c 2 = §. However, one can only combine two 
of the mentioned additional conditions. Condition 9. together with 10. yields 
c 2 = 1 and c 3 = Ci, condition 9. together with ffl~6]) yields c 3 = ±y| while 
condition 10. together with ([TBI) results in c 3 = ci|. One can easily verify that 
all the order 2 conditions 8.-50. are fulfilled with the exception of conditions 
11. and 13.-16. and condition 9. or 10. if (fT6]) is fulfilled in combination with 
only one of them. Therefore, we consider the additional condition (TIB"]) which 
is fulfilled for c\ = 1 and c 2 = c 3 = |. This leads to the SRK scheme RDI1WM 
presented in Table [TJ which is an improved Euler-Maruyama scheme with two 
evaluations of the drift and one of the diffusion coefficients for each step. Thus, 
it may be superior to the widely used Euler-Maruyama scheme, especially in 
practical applications where small noise is inherent to the system. 



4-2 Coefficients for Optimal SRK Schemes of Order (2,2) and (3,2) 

If we consider the order three tree (0^ , o~j 2 , {r, o~j 4 }j 3 ) (see [T0~][l~2"] for details) 
in the case of j\ = j 2 = J3 = J4, then we obtain the corresponding order 
condition 

f3^ T ((A^e)(BWe))(f3^ T e) 2 = \ (22) 
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Table 1 

Coefficients of the optimal SRK scheme RDI1WM with po = 2 and ps = 1. 
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Table 2 

Coefficients of the optimal SRK scheme RDI2WM with po = 2 and ps = 2. 

For the coefficient families f|T0|) -f lT5|) . this order condition is fulfilled if c 3 = 

±v/|- For the tree (0^, {cr^, <Jj 2 , (Jj 3 , Cj 3 }ji) (see [T0lfT2] ) we calculate in the 
case of ji 7^ ji and ji = j'3 the following order three condition 

((3^ T e)((3^ T (B^ef) = 1 (23) 

which is fulfilled if C4 = ±v2. Due to some symmetry in the SRK schemes, 
we obtain always the same SRK schemes regardless what sign we choose for 
c 3 , c 4 and c\. In the following, we choose c 3 = c 4 = and ci = 1. 

Then, the parameter family f lTU]) definitely provides the optimal SRK scheme 
RDI2WM of order pp = p$ = 2 presented in Table [2J 

Next, we calculate SRK schemes of order p D = 3 and ps = 2 for the fa- 
mily (|T3|) in the case (121]) . Again, we try to specify the remaining coefficients 
in the deterministic part of the scheme by additionally considering the order 
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Table 3 

Coefficients of the optimal SRK scheme RDI3WM with p£> = 3 and ps = 2. 
four conditions [3f4] 

a ^(o)(^) e f) = l, (24) 

« T ((A(°)e)(A( )(A(°)e))) = -. (25) 

8 

These conditions are fulfilled if A = | and cs = \. As a result of this, we obtain 
the coefficients of the SRK scheme RDI3WM presented in Table [3j 

However, if we claim for the family ([TBI in the case of (|2"T1) that the order 
four conditions ( jH]) and 

c^e) 3 = 1 (26) 

are fulfilled instead of f )25|) . then we get the coefficients A = 1 and c g = |. As a 
result of this, we obtain the coefficients of the SRK scheme RDI4WM presented 
in Table HJ Here, the deterministic part of scheme RDI4WM coincides with 
the well known Simpson scheme for ODEs [I]. 



5 Numerical example 

In the following, some of the SRK schemes presented in Section H] are applied 
to test equations in order to analyze their order of convergence in comparison 
to some well known schemes. Therefore, the functional u = E(f(X t )) is ap- 
proximated by a Monte Carlo simulation based on the optimal SRK schemes 
RDI1WM of order 1 and RDI3WM and RDI4WM of order 2. The optimal 
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Table 4 

Coefficients of the optimal SRK scheme RDI4WM with pjj = 3 and p,s = 2. 

SRK schemes are compared to the second order SRK scheme PL1WM due 
to Platen [5], the Euler-Maruyama scheme EM of order 1 and the extrapo- 
lated Euler-Maruyama scheme EXEM [5] also attaining order 2. The SRK 
scheme PL1WM is contained in the class of SRK schemes (J3J) with coeffi- 
cients C\ = c 3 = c 4 = 1 in (fit)]) due to case (TA1) . The extrapolated Euler- 
Maruyama approximation is given by 2E(/(Z t V2 )) - E(/(Z*)) based on the 
Euler-Maruyama approximations Z^ 2 and Zj* calculated with step sizes h 
and h/2. The sample average UM,h = jj J2jf=i fO^ti^k)), u k 6 H, of M in- 
dependent simulated realizations of the considered approximation Y t is cal- 
culated in order to estimate the expectation. In the following, we denote by 
A = UM,h — E(/(X t )) the mean error and by a 1 the empirical variance of the 
mean error. Further, we calculate the confidence interval with boundaries a 
and b to the level of 90% for the estimated error fx (see [5fT2] for details). 

As the first example, we consider the non-linear SDE [5|7pTj 

dX t = Ux t + tJx? +T) dt + sfx? + 1 dW t , X = 0, (27) 

on the time interval / = [0, 2] with the solution X t = sinh(t + Wt). Here, we 
choose f{x) = p(arsinh(a;)), where p(z) = z 3 — 6z 2 + 8z is a polynomial. Then 
the expectation of the solution can be calculated as 

E(f(X t )) = t 3 - 3t 2 + It . (28) 

The solution E(f(X t )) is approximated with step sizes 2 _1 , . . . , 2~ 4 and M = 
10 9 simulations are performed in order to determine the systematic error of 
the considered schemes at time t = 2. The results for the applied schemes 
are presented in Table The orders of convergence correspond to the slope 
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Table 5 

Mean errors, empirical variances and confidence intervals for SDE (|27p . 
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of the regression lines plotted in Figure [T] where we get the order 0.58 for 
the EM scheme, order 1.11 for RDI1WM, order 1.80 for EXEM, order 1.81 
for PL1WM, order 1.93 for RDI3WM and order 2.01 for the scheme RDI4WM. 



As a second example, a multi-dimensional SDE with a 2-dimensional driving 
Wiener process is considered: 

d fx}) _(-m ) (xi) dt , ( \xl ±xi \ (wi\ 
[xf) - l-ife -1 + f ) UV + \^x? ±xl + ±x?) {w?) 

(29) 

with initial value Xq = (1, 1) T . This SDE system is of special interest due 
to the fact that it has no commutative noise. Here, we are interested in the 
second moments which depend on both, the drift and the diffusion function 
(see [5] for details). Therefore, we choose f(x) = (x 1 ) 2 and obtain 

E(f(X t )) = exp(-t) . (30) 

We approximate E(f(X t )) at t = 4 by M = 5 • 10 8 simulated trajectories with 
step sizes 2 _0 , . . . , 2~ 3 . The results for the considered schemes are presented 
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Table 6 

Mean errors, empirical variances and confidence intervals for SDE (|29p . 
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Fig. 1. Orders of convergence for SDE ([23) and SDE (j29j) . 

in Tableland Figured) Here, the order of convergence is 0.88 for the Euler- 
Maruyama scheme, 1.53 for RDI1WM, 2.22 for PL1WM, 2.18 for EXEM, 2.24 
for RDI3WM and order 2.28 for the optimal SRK scheme RDI4WM. 

Due to the results in Figure (TJ we can see that for both test equations the 
so called optimal SRK scheme RDI1WM attains much better orders of con- 
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Fig. 2. Computational effort per simulation path versus precision for SDE (p7|) and 
SDE p9j). 

vergence than the well known order one EM scheme. The same holds for the 
optimal SRK schemes RDI3WM and RDI4WM compared to the order two 
schemes EXEM and PL1WM. Clearly, the optimal SRK schemes RDI1WM, 
RDI3WM and RDI4WM need some additional computational effort compared 
to the schemes EM, EXEM and PL1WM, respectively. Therefore, we take the 
number of evaluations of the drift function a and of each diffusion function 
& \ 1 < j — Tn, as well as the number of random variables that have to be 
simulated as a measure for the computational effort. Then we can compare the 
computational effort versus the errors of the analyzed schemes. The results are 
presented in Figure |2j and again, RDI1WM performs much better than the 
scheme EM for both test equations. Further, RDI3WM yields similar results 
like RDI4WM which is for higher precisions slightly better than PL1WM and 
significantly better than the scheme EXEM for the test equation (12"T|) . Con- 
sidering the mult i- dimensional test equation the scheme RDI3WM is 
again close to RDI4WM which performs for higher precisions slightly better 
than EXEM. However both optimal SRK schemes RDI3WM and RDI4WM 
are significantly better than the SRK scheme PL1WM. 



6 Conclusion 



In the present work, a full classification of the coefficients for a class of explicit 
SRK methods of order (1, 1) for s = 1 and order (2, 1) for s = 2 stages as 
well as for the orders (2, 2) and (3, 2) with s = 3 stages is calculated. Based 
on this classification, coefficients for so called optimal SRK schemes are deter- 
mined by considering additional higher order conditions. Optimal coefficients 
for SRK methods of order (2, 1), (2, 2) and (3, 2) are calculated and similarly to 
the deterministic setting [4], better convergence results are expected for these 
schemes in general. Finally, the SRK schemes based on optimal coefficients 
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are applied to some test equations. Here, it turned out that the proposed op- 
timal SRK schemes attain higher orders of convergence than the well known 
schemes under consideration and they also perform very well if the computa- 
tional effort is taken into account. 

For future research, it would be interesting to extend the presented classi- 
fication to diagonal or fully implicit SRK methods. Further, the given classi- 
fication may be applied in order to determine coefficients for SRK methods 
with optimal stability properties. 
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